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ABSTRACT 


In  a  multiple  target  environment  a  radar  signal  processor  often  uses 
weighting  filters  which  are  not  necessarily  matched  to  the  transmitted  wave¬ 
form.  In  this  paper  expressions  for  the  mean- square  ra nge- e stimation  error, 
the  detection  s  ignal- to- noise  (SNR)  and  the  effects  of  sidelobes  are  derived  in 
terms  of  the  impulse  response  of  an  arbitrary  mismatched  filter.  It  is  desired 
to  find  that  impulse  response  which  results  in  the  minimum  range  estimate 
variance  subject  to  preassigned  constraints  on  the  sidelobes  and  the  detection 
SNR.  This  optimization  problem  is  first  formulated  in  state- space  in  which 
the  optimal  control  law  is  sought.  Pontryagin’s  maximum  principle  is  used  to 
obtain  necessary  conditions  for  the  optimum  impulse  response,  from  which  it 
is  possible  to  deduce  the  structure  of  the  optimum  filter.  Certain  mathematical 
details  which  detract  from  the  rigor  of  the  time  domain  formulation  are  resolved 
by  formulating  the  problem  in  the  frequency  domain  and  applying  Hilbert  space 
techniques.  It  is  shown  that  for  the  problem  of  detecting  the  radar  target  and 
estimating  its  range,  the  optimum  filter  is  a  modified  transversal  equalizer. 

If  only  the  detection  function  is  to  be  performed  the  optimum  filter  reduces  to 
the  transversal  equalizer.  This  establishes  the  optimality  of  this  important 
practical  device  as  the  solution  to  the  radar  detection  problem  in  a  multiple 
target  environment.  The  tap  weights  and  spaces  of  the  delay  line  as  well  as 
certain  other  parameters  upon  which  the  solution  depends  can  be  found  by 
solving  a  non-linear  programming  problem.  Numerical  results  are  given  for 
an  interesting  class  of  transmitted  waveforms  which  shows  the  tradeoffs  of  the 
various  filter  parameters. 

Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 
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Introduction 


It  has  long  been  known  that  the  matched  filter  represents  the  "optimum" 
means  of  processing  data  for  obtaining  estimates  of  target  range.  Optimality 
in  this  case  means  that  the  estimates  have  the  smallest  mean- square  error 
possible,  a  re  suit  which  is  valid  only  when  the  si gnal- t o- noi se  ratio  is  large. 

In  many  applications  target  re  solvability  is  a  consideration  almost  as  important 
as  range  accuracy,  and  it  is  then  not  clear  that  the  matched  filter  is  the  best 
receiver  to  use.  In  fact,  for  some  waveforms,  a  sinusoidal  pulse,  for  example, 
very  good  range  accuracy  can  be  obtained,  but  the  resolution  is  poor  because 
the  envelope  of  the  matched  filter  output  signal  has  large  subsidiary  sidelobes. 

One  approach  to  this  problem  is  to  assume  that  the  receiver  is  a  matched 
filter  and  then  try  to  design  the  input  signal  that  will  produce  good  range  estimates 
subject  to  constraints  on  the  sidelobe  structure  of  the  compressed  pulse. 
Algorithms  are  now  available  to  generate  the  solution  to  this  problem  [l],  but 
in  most  cases  the  resulting  waveform  is  quite  complicated,  and  it  is  difficult 
to  build  the  matched  filter  for  it. 

The  other  approach  is  to  use  a  known  signal  which  can  be  transmitted 
easily,  and  to  use  a  mismatched  filter.  This  has  been  done  for  the  linear  FM 
waveform  [2]  in  an  effort  to  reduce  the  sidelobes  of  the  compressed  pulse,  but 
at  the  expense  of  a  loss  in  signal  detectability  and  range  accuracy.  Heretofore, 
no  effort  has  been  made  to  design  a  mismatched  filter  to  simultaneously  minimize 
the  mean- square  range  error  and  loss  of  signal  detectability  subject  to  pre¬ 
assigned  constraints  on  the  sidelobe  structure. 
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In  this  paper  we  assume  that  a  given  pulse  is  received  in  the  presence  of 
additive  white  Gaussian  noise  and  passed  through  a  filter  which  is  not  necessarily 
matched  to  the  input  pulse.  We  assume  that  the  range  is  estimated  by  locating 
the  time  at  which  the  envelope  of  the  filter  output  is  maximum.  The  performance 
of  this  estimation  scheme  has  been  analysed  previously  [3]  with  respect  to 
measuring  the  loss  of  accuracy  and  detectability  due  to  non- optimum  filtering. 

It  is  shown  that  one's  inability  to  build  perfectly  matched  filters  does  result  in 
a  loss  in  performance.  However,  the  advantages  of  mismatching  with  respect 
to  improving  the  multiple  target  resolution  is  not  discussed. 

In  Section  II  we  give  a  brief  derivation  of  the  equations  for  the  mean- square 
range  estimation  error  in  terms  of  an  arbitrary  filter  impulse  response.  The 
effects  of  multiple  targets  are  discussed  and  design  constraints  which  provide 
for  good  resolution  are  developed.  Since  the  processor  is  no  longer  matched 
to  the  transmitted  signal  a  degradation  in  detection  signal- to- noise  must  be 
expected.  (The  detection  signal- to- noise  ratio  for  the  mismatched  filter  is 
derived).  Therefore,  we  define  the  cost  function  as  a  weighted  combination  of 
the  range  estimate  variance  and  the  detection  noise  -  to- signal  ratio.  This 
quantity  is  to  be  minimized  subject  to  the  sidelobe  constraints.  The  weighting 
in  the  cost  function  can  be  varied  to  emphasize  either  the  detection  performance 
or  the  estimation  performance  of  the  filter.  This  weighting  factor  shows  up  as 
a  key  parameter  in  the  structure  of  the  optimal  filter.  We  can  then  associate 
the  detection  and  estimation  functions  of  the  filter  with  its  various  substructures. 

The  filter  design  problem  is  then  formulated  as  an  optimal  control  problem 
in  state  space.  By  following  the  recipe  described  by  Pontryagin's  maximum 
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principle  we  can  find  the  structure  of  the  optimal  filter.  Certain  mathematical 
details  concerned  with  the  stability  of  the  resulting  structure  on  an  infinite 
time  interval  detract  from  the  usefulness  of  the  state-space  technique.  By 
formulating  the  problem  in  the  frequency  domain  and  using  the  insights  obtained 
from  the  optimal  control  solution,  a  mathematically  rigorous  derivation  of  the 
optimum  filter  is  obtained  using  the  projection  theorem  in  a  suitably  defined 
Hilbert  Space.  The  optimum  filter  is  shown  to  be  a  modified  transversal 
equalizer.  If  the  weighting  factor  were  chosen  to  weight  only  the  detection 
signal- to- noise  ratio  in  the  original  cost  functional,  the  optimum  filter  reduces 
to  the  transversal  equalizer.  Therefore,  the  optimality  of  this  structure  for 
target  detection  is  established  rigorously  for  the  first  time  in  the  published 
literature.  For  weighting  factors  which  include  some  measure  of  the  estimation 
performance  it  is  further  established  that  the  transversal  equalizer  is  not  the 
optimal  processor  but  it  is  an  integral  part  of  the  optimum  filter  structure. 

Certain  constants  remain  to  be  determined  to  completely  specify  the 
optimum  filter.  For  example,  the  tap  gains  and  spaces  of  the  transversal 
equalizer  have  to  be  determined  for  a  particular  sidelobe  constraint  function. 

It  is  shown  that  these  constraints  can  be  solved  using  nonlinear  programming 
techniques.  To  show  the  feasibility  of  the  method  and  the  tradeoffs  of  the 
various  filtering  parameters  we  conclude  the  paper  with  a  numerical  evaluation 
of  the  filter  for  an  interesting  class  of  transmitted  waveforms. 
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II  The  Suboptimal  Signal  Processor 


In  this  section  we  shall  derive  the  performance  criteria  to  be  used  in  the 
filter  design  problem.  We  shall  first  concern  ourselves  with  the  variance  of 
the  range  estimate  when  a  mis- matched  filter  is  used  in  the  usual  matched  filter 
or  maximum- likelihood- estimator  configuration.  The  derivation  closely  follows 
that  in  reference  [3]  but  is  included  here  to  make  clear  the  essential  assump¬ 
tions  as  these  become  important  when  the  criterion  is  to  be  used  for  filter 
synthesis.  The  system  analysed  is  shown  in  Figure  L  The  input  consists  of 
a  bandpass  signal  S(t)  plus  white  stationary  bandpass  Gaussian  noise  N(t).  The 
filter,  with  the  transfer  function  Q(w)  =  H(uu-uuc)  +  H(U}-Rjoc),  is  square-law  detected. 
The  time-delay  corresponding  to  the  estimate  of  radar  range  is  then  taken  as 
the  time  at  which  this  envelope  assumes  its  maximum  value.  If  the  filter  were 
matched  to  the  signal,  then  this  processor  is  precisely  the  maximum  likelihood 
estimator  and  has  certain  optimal  properties.  Since  the  filter  is  not  matched  to 
the  signal  the  processor  is  merely  one  specific  suboptimal  method  of  estimating 
radar  range,  (it  is  widely  used  in  practice). 

The  complex  envelope  of  the  input  signal  is 

r  (t)  =  AeJ'9p(t-TQ)  +  n(t)  (1) 

where 

A  =  the  signal  amplitude 

p(t)  =  the  complex  signal  envelope 

9  =  the  random  phase  of  the  received  signal 
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n(t)  =  the  complex  noise  envelope. 

We  let  h(t)  be  the  complex  envelope  of  the  impulse  response  of  the  filter.  Then 
the  signal  at  the  output  of  the  filter  is 

00 

x(t)  =  [  h(r-t)  r(t)dt 
-00 


00  00 

=  A€J0  J  h(T-t)p(t-To)dt  +  J  h(T-t)n(t)dt 
-  00  -  00 


We  define  the  quantities 


(2) 


c p(r)  =  J*  h(r-t)p(t)dt  (3a) 

_  00 

00 

T|(t)  =  J  h(T-t)n(t)dt  (3b) 

-  00 

and  then  write  for  the  filter  output 

x(t)  =  Ae^0CD(T-To)  +  T|(T)  (4) 

It  is  at  this  point  that  our  analysis  first  begins  to  digress  from  that  in  reference 
[3].  It  is  extremely  convenient  for  the  work  which  follows  to  restrict  the  class 
of  waveforms  and  filters  to  be  those  for  which  the  filtered  signal  (which  we  may 
sometimes  refer  to  as  the  compressed  pulse)  is  real.  This  may  appear  to  be  a 
very  restrictive  assumption,  but  in  fact  it  is  not.  It  turns  out  that  after  we  have 
solved  for  the  optimum  mis-matched  filter  this  assumption  is  equivalent  to  re¬ 
quiring  that  the  autocorrelation  functions  of  the  input  waveforms  be  real.  This 
in  turn  will  be  valid  for  signals  having  symmetric  spectral  densities,  which  is 
the  case  of  most  practical  interest.  For  example,  the  linear  FM  waveform 
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falls  within  this  class.  We  could  have  refrained  from  making  this  assumption 
until  a  later  point  in  the  paper,  but  we  felt  that  since  we  were  going  to  make 
the  assumption  anyway,  we  might  as  well  do  it  earlier  in  the  paper  since  it 
makes  the  manipulations  at  this  point  quite  a  bit  easier.  Therefore  the  output 
of  the  envelope  detector  can  be  written  as 


z(T)  =  |  A£J9CD(T-To)  +  T\(t)\ 


2 


=  A2cp2(T-To)  +  2Aco(r-To)Re[.e^'n(T)]  +  |n(r)|“  ^ 

The  time-delay  cor  re  spondi  ng  to  the  target  range  is  taken  as  the  time  at  which 

z (t)  achieves  its  maximum  value.  In  other  words  the  estimate  of  r  is  the  number 

o 

T  where 

v  max  .  ,  /c . 

z(t)  =  T  z(t)  (5) 

In  the  matched  filter  case,  h(t)  =  p* (- 1 ),  and  in  the  absence  of  noise,  it  follows 

from  (3a)  that 


z(r)  =  [  J  p(t)p*(t-r+ro)dtr  (6) 

_  00 

which  is  the  square  of  the  autocorrelation  function  of  the  transmitted  pulse.  This 
quantity  is  indeed  maximized  at  the  true  parameter  value  T  .  In  the  mis-matched 
case  and  no  noise 


z(t)  =  A2[cd(t-to)]2  (7) 

The  function  c p(T)  will  achieve  its  maximum  value  at  a  point  which,  in  general, 

may  be  non-zero.  In  this  case  it  becomes  clear  that  we  should  choose  as  our 

estimate  of  T  the  number  f,  -T.  where 
o  b  b 


7 


Z<V  = 


CD(Tb)  = 


max 


max 


z(T) 


cd(t  ) 


(8a) 

(8b) 


It  is  possible  to  develope  criteria  describing  the  detection,  estimation  and  re¬ 
solution  capabilities  of  the  mis-matched  filter  in  terms  of  an  arbitrary  bias 
T^.  However  it  can  be  shown  that  enters  the  performance  equations  only 
through  the  term  h(t-T^).  Consequently  the  performance  of  the  optimum  biased 
filter,  h^(t),  can  always  be  achieved  by  using  an  unbiased  filter  h^(t)  by  setting 
h^(t)  =  h^(t-T^),  that  is,  by  inserting  units  of  delay  in  cascade  with  h^(t). 
Therefore,  with  no  loss  in  generality  we  can  limit  our  attention  to  those  filters 
for  which  the  bias  is  zero.  Hence  we  need  consider  only  filters  for  which 
Co'(0)  =  0,  (the  prime  denotes  the  derivative  with  respect  to  the  argument). 

The  remainder  of  the  analysis  proceeds  once  again  along  the  lines  of  that 
described  in  [3],  We  want  to  point  out,  however,  that  the  motivation  in  that 
paper  was  to  determine  the  loss  in  accuracy  as  a  result  of  using  an  imperfectly 
constructed  matched  filter.  No  mention  is  made  there  of  the  idea  of  synthesizing 
a  mis-matched  filter  to  give  good  range  estimates  and  low  sidelobes  simul- 
ta  neously. 

We  shall  now  return  to  the  immediate  problem  of  determining  the  equations 
describing  the  range  variance  of  the  ad  hoc  estimator.  First  we  write  the 
output  of  the  envelope  detector,  z (r),  as  a  Taylor  series  in  T  about  the  point  r  , 


z(r)  =  z(r  )  +  z'(r  )(r-r  )  +  y  z"  (r  )(t-t  )2  +  .  .  . 
o  o  o  Z  o  o 


(9) 


For  large  signal- to- noise  ratio  the  estimate  T  will  be  ’’close  to”  the  true  value 


Tq  and  we  can  expect  that 
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(10) 


z(f)  =  z(r  )  +  z ' (t  )(f-r  )  +  i  z"  (r  )(f-r  )2 

'  o  o  o  2  o  o 

But  f  is  the  value  which  maximizes  z  (t),  that  is,  z  7  (f )  =  0.  Therefore  it 
follows  that  the  estimation  error  is 

T-T  =  -z'(T  )/z//(r  )  (11) 

o  o  o 

We  eventually  want  an  expression  for  the  error  which  is  to  first  order  in  the 
noise.  This  can  be  obtained  by  noting  that  for  large  SNR  the  second  order  noise 
term  in  (4)  can  be  neglected  so  that 

z(r)  =  A^2(t-tq)  +  2Acc  (T-T^)Re[c^Ari(T)]  (12) 

and  therefore  we  can  write  for  z  '  (r  )  the  equation 

z ' (t)  =  2Ac c'(T-to)  |acp(t-T(i)  +  Re[V9ri(T )]  )  (13) 

+  2Acd (r-ro)Re[e^ri '  (r)] 

Since  all  filters  under  considera tion  are  unbiased,  cd^O)  =  0,  and  therefore 

z'(to)  =  2Acp(0)Re[-pj9n'(r)]  (14) 

This  renders  the  numerator  in  (11)  to  be  first  order  in  the  noise  whence  it 
suffices  to  express  the  denominator  to  zero  order  in  the  noise.  Therefore  we 
have  from  (13) 

z"(ro)  =  2AZcd(0>£/,(0)  (15) 

and  the  estimation  error  is 

T-r  =  -  Re[ej  V(to)]/Acc"  (0)  (16) 
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From  (3b)  it  follows  that 


r1/(ro)  =  J  h'(To-t)n(t)dt 


which  is  a  zero  mean,  Gaussian  random  variable.  For  the  complex  noise 
envelope  one  has^ 


n(t)n(s )  =  0 


IvfFJn^TsT -  2NQ5(t-s) 


whence  it  follows  that 


t1,(To)ti,(to)  =  0 


ti,(to)ti/*(to)  =  2Nq  J  |h'(ro-t)|  dt 

_  oo 

Then  if  we  define  the  complex  variable 


(17a) 

(17b) 


C  =  C  +  jC  =  Jeri'(T  ) 
x  J  y  !  x  o 


it  follows  from  (17a)  that 
”2 


c  =  C  -C  +2jCC  =0 

x  y  J  x  y 


so  that  C  C  =0  and  C  =  C  Using  these  facts  and  (17b)  it  follows  that 
x  y  x  y  ° 


CC*  =  C  +C  =2C  =2Nf  1  h  7  (- 1  )  I  dt 
x  y  x  oJ 

J  m 


the  bar  over  a  quantity  denotes  the  ensemble  average. 

i 

'N  is  the  (si  ngl  e- si  ded )  noise  spectral  density. 
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Therefore  the  estimation  error,  Equation  (16),  has  zero  mean  and  variance 


(f-T 


N 

o 


.  J  |h'(-t)|2dt/fc"(0)] 

_  00 


From  (3a)  we  find  that 


(IB) 


00 

co"  (0)  =  t  h"  (- t)p(t)dt 

_  oo 


(19a) 


=  -  j  h'(-t)p'(t)dt  (19b) 

_  00 

where  (19b)  follows  from  an  integration  by  parts.  Combining  (18)  and  (19)  we 
conclude  that  for  large  SNR  the  mis-matched  filter  leads  to  unbiased  estimates 
of  the  true  time  delay  r  .  The  variance  of  the  estimate  is  given  by 


^  2 
mmf 


J  ih'(-t)i2dt 


00 

J  h,(-t)p/(t)dt 
-00 


2 


(20) 


and  this  result  is  vralid  provided  the  signal  and  filter  result  in  a  compressed 
pulse,  cp  (t  ),  which  is  real  and  has  a  maximum  value  at  the  origin,  namely 
cp  '(0)=  0  or 


h  '(- t)p(t)dt  =  0 


(21) 


If  the  filter  is  matched  to  the  signal,  h(t)  =  then  (21)  is  satisfied  and 


(20)  becomes 

00 

°  mf  =  V  I  |p'<t»|2dt  <22> 

-  00 

which  is  the  classical  result.  It  is  well-known  [4]  that  the  matched  filter  pro- 

2  >2 

duces  the  smallest  mean- squared  ranee  error  so  that  a  r  =  a  r.  Then 

^  °  mmf  mf 


2  J  |h'(-t)|2dtj  |p'(t)|2dt 

Q  „  00  _  00 

mmf 

2  00 
mf  [  J  h'  (-t)p'(t)  dt]2 


(23) 


measures  the  increase  in  the  e  stimation er r or  as  a  result  of  using  a  mismatched 
filter . 

So  far  we  have  considered  only  the  detrimental  effects  of  mismatching  in 
relation  to  the  degradation  in  range  performance  as  expressed  in  Equation  (23). 
The  major  reason  for  using  deliberate  mismatching  is  to  achieve  better  per¬ 
formance  in  the  multiple- tar  get  environment.  In  a  typical  application  the 
waveform  CP(t)  will  have,  in  addition  to  a  maximum  at  the  origin,  subsidiary 
peaks,  called  sidelobes,  at  other  values  of  T.  If  these  sidelobes  are  not 
significantly  smaller  than  the  mainlobe,  then  in  the  multiple-tar  get  environment 
it  is  possible  that  targets  in  adjacent  range  cells  with  smaller  radar  cross- 
sections  will  be  masked  by  the  sidelobes  of  a  large  radar  cross-section  target. 

It  is  desirable,  therefore,  to  design  the  filter  to  make  these  sidelobes  small 
with  respect  to  the  magnitude  of  the  central  lobe  at  T  .  This  can  be  accomplished 
by  requiring  that  the  constraint 


We  can  assume  that  cp(  0)  >  0  since  h(t)  can  always  be  replaced  by  -h(t)  without 
changing  the  estimation  performance  or  the  zero-bias  constraint. 
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-€(t)cp(0)  =cp (r)  =  €(r)cp(0) 


(24) 


be  satisfied,  where  0  =  6(r)  =  1  represents  the  sidelobe  constraint  function  and 
is  chosen  to  combat  the  particular  clutter  distribution  under  consideration.  In 
many  practical  problem  nothing  is  known  concerning  the  clutter  distribution  to 
be  encountered  and  a  useful  constraint  function  is  to  set  £(t)  =  1  for  |r|<  6, 

€(t)  =  (1-y)  for  |t|>  6  where  *  is  the  width  of  a  range  resolution  cell  and  y  may¬ 
be  from  .  9  to  .  99  in  extreme  cases.  Notice  that  (21)  in  conjunction  with  (24) 
guarantees  that  the  extremum  of  cp(T)  at  the  origin  is  a  global  maximum.  Hence¬ 
forth  we  shall  require  (with  no  loss  in  generality)  that 

cp ( 0 )  >  0  (25a) 

-CP  (0)  >0  (25b) 

Using  (3)  in  (24),  the  constraint  on  the  structure  of  the  compressed  pulse  will 
be  satisfied  provided 

00  00 

|J  h(-t)p(t+r)dt|  =  €(t)  J  h(-t)p(t)dt  (26) 

-  00  _  00 

In  practice  h(t)  and  p(t)  will  be  bandlimited  functions,  in  some  sense  and  there¬ 
fore  "slowly  varying"  functions  of  T.  Therefore,  the  continuum  of  constraints 
can  be  replaced  by  the  finite  number  of  constraints 

< 


h(-  t)p(t+r  .  )dt  |  =  G.  J  °°h(-  t)p(t )dt  j  =  l,  2 

J  J  on 


.  .  .  ,  n 


(27) 


where  nand  {r.}  are  chosen  to  provide  a  dense  sampling  of  the  compressed  pulse. 

Jj=i 

The  problem  for  which  this  approximation  is  not  made  has  been  studied  recently 
in  reference  [5].  It  is  shown  there  that  the  mathematical  framework  needed 
to  handle  the  continuum  of  constraints  is  so  complicated  that  one  not  only  loses 
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touch  with  the  physical  problem  at  hand  but  also  finds  it  impossible  to  interpret 

the  solution  once  it  is  obtained.  We  shall  see  that  the  discretized  problem  can 

readily  be  solved  and  interpreted  and  we  can  then  show  that  as  the  sampling 

becomes  dense,  the  sequence  of  optimum  solutions  converges  with  respect  to 

an  appropriately  defined  Hilbert  Space  norm. 

In  order  to  complete  the  specification  of  the  filter  design  problem  we  must 

take  into  account  the  loss  of  signal  detectability  which  must  occur  as  a  result 

of  the  filter- signal  mismatch.  That  such  a  loss  occurs  follows  from  the  well- 

known  fact  that  the  matched  filter  results  in  the  largest  detection  signal- to- noise 

ratio  (SNR).  In  order  to  develop  a  criterion  which  measures  this  loss  we  shall 

adopt  a  detection  rule  which  is  used  in  the  matched  filter  processor.  That  is, 

a  radar  target  is  said  to  be  encountered  if 

z (T )  >  \  (28) 

and  that  no  target  is  present  otherwise.  In  this  case  T  is  the  estimate  obtained 

by  solving  (5)  using  the  processor  described  above,  and  X  is  a  suitably  chosen 

2 

threshold  level.  Since  cp  (t-T  ),  Equation  (3),  represents  the  instantaneous 
signal  power  at  the  output  of  the  filter,  and  since 

ZN,  J  >(-t),2dt  (29) 

_  00 

represents  the  average  noise  power,  the  detection  noise  -  to- signal  ratio  is 
measured  by  the  quantity 

00  _  _ 

2No  !M-t)|2dt/cp2(f-To)  (30) 

A  well- designed  receiver  will  provide  good  estimates  of  T  so  that  3=- T  ~  0  and 

o  o 

the  detection  performance  of  the  mismatched  filter  is  described  by  the  noise- 
to-  signal  ratio 


14 


J  l^(-t)|2  dt 


p  =  2N 
mmf  ° 


[  h(-t)p(t)dt] 
for  the  matched  filter  case,  h(-t)  =  p*(t),  and  we  have 


(31) 


Pmf  =  2No7  J*  IP(t),2dt 

_  CO 

The  ratio 


(32) 


P 


mmf 


Pmf 


| h (-  t)  |2  dt  lp(t)pdt 
—  00  —  00  > 

-  =  1 

o  00  2 

[  J*  h(-t)p(t)dt] 

_  00 


(33) 


represents  the  increase  in  the  detection  noise- to- signal  ratio  as  a  result  of 
using  a  mismatched  filter.  Since  the  signal  p(t)  is  specified  we  can  define 


J  |P(t)  (2  dt  =  E  (34a) 

-  00 

00 

I  lp'(t)|2dt  =  B2  (34b) 

~°°  z 

E  represents  the  signal  energy  which  B  /E  represents  its  mean  square  bandwidth. 

On  the  basis  of  the  preceding  analysis  we  can  combine  Equations  (21),  (23), 
(27),  (33),  and  (34)  to  formulate  two  important  filter  design  problems:  From  the 
class  of  admissible  filters,  we  want  to  find  that  filter,  which,  for  a  given  signal 
p(t),  satisfies  the  multiple  target  resolution  constraints,  Equation  (27). 


J  h(-t)p(t+T.)dt  |  =  €  J  h(-t)p(t)dt  j  =  l.  2,...,n  (35) 

no  ^  ^  m 
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and  satisfies  the  zero  bias  constraint,  Equation  (21) 


Joo 

h '  (t)p(t)dt  =  0 


(36) 


and  extremizes  the  cost  functionals  for  the  following  problems: 

Problem  1:  minimizes  the  normalized  detection  noise- to- signal  ratio, 
00 

Equation  (33)  |h(- t)  |2dt 

E  •  — - 


[  J*  h(-t)p(t)dt]  (37) 

-  00 

Problem  2:  minimizes  the  normalized  increase  in  the  mean- square  range 
estimation  error  Equation  (23) 


B2  • 


[J  h'(-t)p'(t)dt]‘ 


(38) 


subject  to  a  constraint  on  the  allowable  increase  in  detection  noise- to- signal 


Equation  (33) 


E 


J  N-t)|2dt 
-  00 
00 

[J  h(-t)p(t)dt]‘ 


=  e  >  i 


(39) 


It  is  advantageous  to  combine  these  two  problems  into  a  one- parameter 
family  of  optimization  problems  which  minimizes  a  linear  combination  of  the 
detection  and  estimation  performance  indices.  Therefore,  we  combine  Equatio 
(37),  and  (38)  and  (39)  and  minimize 

rih'Mifd.  . 2 


ns 


a 


B 


J  00  |h(-t)Tdt 


E 


uu 

[J  h'(-t)p,(t)dt]' 


+  (1-a) - — - for  o  =  a  =  1 


[  J  h(-t)p(t)dt]" 


(40) 
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subject  to  the  zero  bias  constraint,  Equation  (36) 


00 

J  h  '  (-t)p(t)dt  =  0 
_  00 


(41) 


and  sidelobe  constraints,  Equation  (35) 

00  00  00 
J  h(-t)p(t+r.)dt  =  J  h(-t)p(t)dt=  €.J  h(- t)p(t+r .  )dt 

J  _  m  J  _  OO  J_00 


j=lt  2,  .  .  .  ,  n  (42) 

By  setting  a  =  0  we  could  obtain  the  filter  which  maximizes  the  detection 
signal- to- noise  ratio  subject  to  sidelobe  constraints.  This  might  be  a  useful 
solution  to  obtain  for  some  communications  problems  in  which  the  signal  must 
"stand  out"  of  the  noise.  For  a  =  1  the  solution  would  yield  a  filter  which  gives 
the  best  range  accuracy  subject  to  sidelobe  constraints.  It  would  be  tempting 
to  use  this  filter  for  pulse  position  modulation  communications  problems  and 
for  the  radar  tracking  problem.  However,  one  intuitively  would  expect  the 
performance  to  depend  on  the  ability  to  detect  the  target,  since  if  the  signal 
does  not  "stand  out"  of  the  noise,  it  would  not  be  possible  to  measure  the 
location  of  the  peak  of  the  compressed  pulse.  In  fact  we  shall  show  that  the 
case  a  =  1  leads  to  a  mathematically  ill-defined  problem.  In  addition,  the 
formula  describing  the  estimation  performance  of  the  filter  implicitly  assumed 
a  large  signal- to- noise  ratio.  This  assumption  can  not  be  guaranteed  unless 
a  <  1. 

For  the  cases  in  which  0  <  a  <  1,  we  are  merely  putting  different  weights 
on  the  detection  and  estimation  performance  characteristics  of  the  filter.  By 
solving  the  problem  for  several  values  of  a  we  can  choose  the  filter  which 
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gives  an  acceptable  detection  performance  and  good  estimation  performance. 

It  can  be  shown  that  that  these  solutions  correspond  to  those  of 
Problem  2  for  appropriate  values  of  f3.  Therefore,  the  problem  which  is  de¬ 
scribed  by  Equations  (40),  (41),  and  (42)  describes  a  wide  class  of  problems 
which  arise  in  radar  and  communications.  It  is  to  this  latter  problem  formula¬ 
tion  that  we  shall  direct  our  attention  in  the  following  sections.  In  the  next 
section  we  shall  reformulate  this  problem  in  state- space  and  by  using  optimal 
control  techniques  we  shall  obtain  necessary  conditions  for  the  optimum  filter. 
From  these  conditions  we  can  deduce  the  physical  configuration  of  the  optimum 
filter.  However,  certain  mathematical  difficulties  arise  in  the  time  domain 
formulation  which  are  difficult  to  resolve.  By  transforming  the  problem  to 
the  frequency  domain  and  using  the  insights  obtained  from  the  structure  of  the 
time  domain  solution,  it  is  possible  to  use  functional  analysis  to  obtain  the 
Fourier  Transform  of  the  unique  optimum  filter.  This  solution  specifices  the 
optimum  filter  except  for  certain  unknown  gains.  The  cost  functional,  the 
zero  bias  constraint  and  the  sidelobe  constraints  can  be  evaluated  in  terms  of 
these  unknown  constants  which  in  turn  leads  to  a  nonlinear  programming 
problem.  Efficient  computational  algorithms  which  have  been  developed  for 
this  type  of  problem  are  then  used  to  generate  numerical  results  to  illustrate 
the  utility  of  this  approach  to  filter  design. 
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Ill,  State-Space  Solution  of  the  Design  Problem 


In  this  section  we  shall  formulate  and  solve  the  mismatched  filter  design 
problem  using  state- va riable  techniques.  Most  of  the  section  involves  straight¬ 
forward  details  which  will  provide  no  new  insights  to  anyone  familiar  with 
optimal  control  theory.  It  is  included  at  this  point  in  the  paper  because  this 
was  the  order  in  which  we  performed  our  research.  It  shows  that  by  following 
the  recipe  prescribed  by  the  minimum  principle  of  Pontryagin  one  can  formally 
obtain  the  structure  of  the  optimum  filter.  With  the  solution  before  us,  we 
were  then  able  to  obtain  a  much  more  rigorous  derivation  of  the  result  which 
uses  only  the  projection  theorem  of  Hilbert  Space  Theory.  The  latter  solution 
is  presented  in  Section  IV.  The  point  is  that  Section  III  can  be  left  out  by  the 
casual  reader  without  significant  loss  in  continuity.  It  is  included  for  those 
who  may  have  an  interest  in  the  methodology  which  must  be  followed  to  use 
optimal  control  theory  to  solve  optimization  problems.  Some  limitations  of 
the  method  are  pointed  out. 

In  order  to  guarantee  that  the  optimal  control  approach  will  result  in  a 
realizable  impulse  response,  we  shall  replace  the  function  h(-t)  which  appears 
in  all  of  our  equations  by  h(T-t)  where  T  denotes  some  suitably  large  processing 
time.  We  shall  assume  that  h(t)  and  p(t)  are  real.  In  addition,  we  assume 
that  p(t)  and  p(t-fT.),  j  =  ;,  .  .  .  ,  n  are  zero  outside  of  the  interval  [0,  T]. 

We  now  proceed  to  the  state- space  formulation  by  defining  the  control 
function  to  be  the  first  derivative  of  the  impulse  response  and  let  the  first 
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component  of  the  state  vector  be  the  impulse  response*,  i.  e. 


u(t)  =  h(T-t)  (46a) 

Xj(t)  =  h(T-t)  (46b) 

Then  the  first  state  equation  is 

x^t)  =  -  u(t)  (47) 

In  order  that  H  (T-t)  be  well  defined  for  all  te[0,  T]  we  require  that  x^(0)  =  x^(T)  = 
0.  In  addition  we  define  the  state  variables 

i^(t)  =  p(t)K  (T-t)  (48a) 

x3(t)  =  p(t)K(T-t)  (48b) 

i4(t)  =  h2(T-t)  (48c) 

i5(t)  =  li2(T-t)  (48d) 

yQ(t)  =  p(t)h( T-t)  (48e  ) 

Yj(t)  =  p(t+r.)h(T-t)  j  =  l,  2,...,n  (48f) 

Each  of  these  variables  is  to  have  zero  initial  conditions,  whence  it  follows  that 

T 

x2  (T)  =  J  p(t)K(T-t)dt  (49a) 

o 

T 

x3(T)=J  p(t)H(T-t)dt  (49b) 

o 

T  2 

x4(T)=J  [h(T-t)]  dt  (49c) 

o 

T 

x5(T)  =  J  [K(T-t)]2dt  (49d) 

o 


*  The  dot  denotes  differentiation  with  respect  to  time. 
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(49e) 


T 

y  (T)  =  J  p(t)h( T- t)dt 

o 


T 

y  (T)=J  p(t+T  )h(T-t)dt  (49f) 

J  o  ** 

j  =  l,  2,  .  .  .  ,  n 


Notice  that  these  are  all  of  the  quantities  needed  to  specify  the  filter  design 
problem.  By  substituting  Equation  (46)  into  Equation  (48)  we  obtain  the  state 


equations 

x^t)  =  -  u  (t )  x^O)  =  0  XjjT)  =  0  <50a) 

x2(t)  =  p(t)u(t)  x2(0)  =  0  (50b) 

Xj(t)  =  p(t)u(t)  x3(0)  =  0  (50c) 

x4(t)  =  x^(t)  x4(0)  =  0  (50d) 

x5(t)  =  u2(t)  X5  (0)  =  0  (50e ) 

yQ(t)  =  p(t)x1(t)  yQ(0)  =  0  (50f ) 

y.(t)  =  pft+r.lx^t)  y^(0)  =  0  j  =  l,  2,  .  .  .  ,  n  (50g) 

The  resolution  constraints,  Equation  (42),  require  that 


-€  y  (T)  =  y.(T)  =  €.y  (T)  j=l,  2 . a  (51) 

Jo  J  jo 

and  the  zero  bias  constraint,  Equation  (41),  requires  that  x^(T)  =  0.  The  cost 
functional  to  be  minimized,  from  Equation  (40),  can  be  written  as 


x5  (T) 

P[x(T)]  =  a  .  -^7— 

x2 


+  (1-a)  . 


x4(T) 

*7™ 


(52) 


Since  all  of  the  equations  are  homogeneous  in  h(*  ),  we  can  normalize  either 
x?(T)  or  y  (T).  Therefore,  the  filter  design  problem  is  equivalent  to  finding 

L>  O 

the  control  function  u(t),  te[o,  T]  which  minimizes  the  cost  functional 


x5m  r2 

P[x(T)]  =  a - 2—  •  +  (1-a)  x  (T)  (53) 

E  4 

subject  to  the  differential  equation  constraints 

Xj(t)  =  -u(t)  x1(0)  =  0  Xj(T)  =  0  (54a) 

x2(t)  =  p(t)u(t)  x2(°)  =  0  (54b) 

x3(t)  =  p(t)u(t)  x3(0)  =  0  x3(T)  =  0  (54c) 

x4(t)  =  X^t)  x4(0)  =  0  (54d) 

k5(t)=u2(t)  x5(0)  =  0  (54e ) 

yQ(t)  =  p(t)Xl(t)  yQ(0)  =  0  yo(T)  =  1  (54f) 

Yj(t)  =  P(t+T^)x^(t)  y  j  ( 0 )  =  0  j=l,  2,...,n  (54g) 


subject  to  the  resolution  constraints 


=  Yj ( T)  =  6.  j=l,  2,  ....  n 


(55) 


Although  the  dynamics  of  the  system  equations  are  relatively  simple,  the 
inequality  end-point  constraints,  Equation  (55),  makes  this  an  interesting 
optimal  control  problem.  We  shall  deal  with  these  constraints  using  the  method 
suggested  in  [6]  in  which  we  imbed  the  original  problem  into  a  subclass  of  pro¬ 
blems  in  which  the  terminal  conditions  are  fixed.  Therefore,  we  shal]  solve  a 
fixed  end-point  problem  by  requiring  that 

y.(T)=y.f  j  =  l,  2,  .  .  .  ,  n  (56) 
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(57) 


where  the  n- vector  belongs  to  a  set  S  c  En  where 


s  =  (x1:  j  =  l.  2,...,n} 

For  each  vector  we  have  a  fixed  terminal  time  optimal  control  problem 

for  which  the  necessary  conditions  for  optimality  can  be  derived  using  the  min- 
principle  [7].  Consider  the  general  control  system  described  by 

z(t)  =  fjz(t),  u(t),  t]  0  =  t  =  T  (58a) 


imum 


z(0)  =  z 
g.[z(T)]  =  0 


j  =  l,  2,  .  .  .  ,  r 


(58b) 

(58c) 


in  which  a  minimum  of  P[z(T)]  is  sought,  we  first  form  the  Hamiltonian  function 

N 

H[z(t),  _X(t),  u(t),  t]  =  .  |  l  \.(t)f.[z(t),u(t),  t]  (59) 

where  the  costate  variable  \(t)  satisfies  the  differential  equation 


X*  .(t) 

with  e nd- conditions 


3H 


?Zj(t) 


j  =  1.2,  ...  ,  n 


(60a) 


Xj(T) 


3P[z(T)] 

SZj(T) 


+  .Z.a. 

1=  1  1 


9gJz(T)] 

AZj(T) 


(60b) 


The  constants  are  unknown  multipliers  which  are  to  be  chosen  so  that  the 
prescribed  end- conditions,  Equation  (58c),  will  be  satisfied. 

Applied  to  the  problem  at  hand  the  Hamiltonian  is, 


H  =  \1(t)u(t)  +  \2(t)p(t)u(t)  +  \3(t)p(t)u(t) 

+  X4(t)x1Z(t)  +  X5(t)u2(t)  +  UQ(t)p(t)x1(t) 


1=1  uk(t)P(t+Tk)x1(t) 


(61) 
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Therefore,  the  costate  variables  satisfy  the  equations 

n 


Vj(t)  =  -  2X4(t)x1(t)  -  Uo(t)p(t)  -  E  uk(t)p(t+rk) 

k=l 

X.(t)  =  0  j=2,  .  .  .  ,  5 

Ctj  (t)  =  0  j=0,  1,  . . .  ,  n 

The  end-point  constraint  functions  are 

gl[z(T)]  =  Xj(T) 

g2[z(T)]  =  x3(T) 


g3[z(T)]  =  yo(T)  -  1 


g3  +j  — 

and  the  cost  functional  is 


,[z(T)]  =  y.(T)  -  y.1  j=l,2,...,n 


P[z(T)]  =  a 


x5(T) 

x92(T) 


B 


E  +(l-a)x4(T) 


Applied  to  Equation  (60b),  the  terminal  values  of  the  costate 


Xj(T)  =  Oj 

X2(T)  =  -2a 


x5(T)  b2 


x2'(T) 


X3(T)  =a2 


X4(T)  =  (1-0.) 


X5(T)  =o  .—4—  •  % 

x2  (T) 


U  (T)  =  a, 
o  3 


aj(T)  =  a3+^  j=l,  2,  ....  n 


(62a ) 
(62b) 
(62c) 

(63a) 

(63b) 

(63c) 

(63d) 

(64) 

variables  are 

(65a) 

(65b) 

(65c) 

(65d) 

(65e ) 
(65f) 
(65g) 
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From  Equation  (62b),  (62c),  it  is  clear  that  all  costate  variables  except  ) 
are  constants.  Therefore,  the  Hamiltonian  can  be  written  as 


xc(  T)  r2 

H  =  [-X  (t)  -  2tt  .  -4^—  .  p(t)  +  a?P(t)]  u(t) 

1  x23(T)  E  2 

+ — ^ u2(t)  +  (l-a)x^2(t)  +a3p(t)x1(t) 
x2  ^T) 

+  kE=lCt3+kp(t+Tk)  Xl(t)  (66) 


The  minimum  principle  states  that  the  optimal  control  must  minimize  H  at  each 
t*[0,  T].  The  minimization  is  performed  over  the  class  of  admissible  controls, 
which  in  turn,  is  directly  related  to  the  class  of  admissible  filters.  If  any 
nbuildabi  li  ty n  constraints  are  to  be  incorporated  into  the  synthesis,  it  is  at 
this  point  that  they  are  to  be  taken  into  account.  For  the  purpose  of  this  in¬ 
vestigation,  however,  we  shall  assume  that  there  are  no  constraints  on  the 
filter's  structure.  Then  Equation  (66)  can  be  minimized  over  the  variable 
u(t)  for  each  t~[0  T]  by  minimizing 


*c(T) 

H[u(t)]  =  [  -  X ,(t)  -  2a  p - 

x2  <T> 

,2 


B 


V 


-p(t)  +  a9p(t)]  u(t) 


a 

"Z 


(T) 


(67) 
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Recall  the  fact  that  the  weighting  constant  &e[0,  l].  For  0<a  =  1,  the  function 
H  has  a  well-defined  minimum  which  is  achieved  by  the  control 


x22m 


F  r  xc(T) 

.  — 2  •  [  ,(t)  +  2a  — ^ —  .  p(t)  -  a^p(t)] 

-2  <T) 


(68) 


B 


However  if  a  =  0,  then  the  control  is  singular  [8]  provided 


XJ(t)  -  a2p(t)  h  o  for  all  t€[Q,  T] 


(69) 


Otherwise  the  optimum  control  switches  from  ±°°  which  is  a  meaningless  solution. 
Information  regarding  the  nature  of  the  optimum  singluar  control  can  be  obtained 
by  differentiating  Equation  (69)  and  relating  the  result  to  Equation  (62a).  We 
obtain 


a2p(t)  =  X*  ^t) 

n 

=  -2xj*(t)  -  UQp(t)  -  kE1ukp(t+Tk) 


(70) 


where  we  have  set  =  01^^,  k  =  0,  1,  .  .  .  ,  n.  But  by  definition  x^(t)  =  h(T-t), 
therefore,  the  optimum  impulse  response  is  of  the  form 

n 

h*(T-t)  =  X  xp  (t )  +  X2p(t)  +kE^kp(t+rk)  (71) 


where  \ ,  U-,  .  .  .  , U  are  arbitrary  constants.  These  are  to  be  chosen  to 
satisfy  the  sidelobe  and  zero-bias  constraints  and  to  further  minimize  the  cost 
function.  We  shall  discuss  this  idea  at  length  at  a  later  time.  The  point  we 
want  to  make  here  is  that  the  optimal  control  methodology  has  led  us  to  the 
form  of  the  optimum  filter,  at  least  for  the  case  a  =  0.  But  setting  a  =  0  in 
the  cost  functional,  Equation  (53),  means  that  no  consideration  was  being  given 
to  the  range  estimate  variance;  only  the  detection  noise- to- signal  ratio  was 
being  minimized.  Therefore,  the  singular  control  leads  to  the  structure  of  the 
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filter  which  maximizes  the  detection  signal- to- noise  ratio  subject  to  fixed  con¬ 
straints  on  the  sidelobe  structure.  Equation  (71),  therefore,  proves  that  the 
transversal  equalizer  is  the  optimum  filter  for  maximizing  the  detection  SNR 
in  a  multiple  target  environment. 

Next  we  return  to  the  case  0<CC  =  1.  For  this  case  the  optimal  control  is  given 
by  Equation  (68).  However,  since  x^(t)  =  -u(t),  then  Equation  (68)  can  be 
written  as 


(t) 


E  X22<T) 


\(t)  - 


x5  ( T) 

3^rry 


p(t)  + 


x /  (T) 

-TT-PW 


(72) 


We  define  a  new  function  according  to  the  equation 

*  *5<T) 

xi’(t)  =  e(t)  -  ihjTj  p(t)  (73) 

which,  when  substituted  in  Equation  (72)  leads  to 

V  X?Z(T)  F  X22(T)  . 

^(t)  "  '  — Za — *  ^i<t)  +  ~2  Za — a2p(t)  (74) 

B  B 


Since  the  derivatives  of  X^(«  )  and  p(*  )  have  been  assumed  to  exist,  then 


E  x22(T) 


E  X2  <T> 


® (t) - z  ’  ~2a -  ‘  X  I(t)  +_Z  ““Za —  a2p(t) 


B 


B 


and  then  using  Equation  (62a)  for  X  ^(*)  ,  we  have 


x  2(T)  n 

e(t)  =  -  — 2  — [-2(l-a)Xl  (t)-uop(t)-  k|j  ukP(t+Tk)] 

B 

F  X?2(T> 

+  rz  “2a -  a2p(t) 

JD 


(75) 


(76) 


Then  utilizing  Equation  (73)  for  x^(t)  we  have  the  second  order  differential 
equation  for  e(*  ), 
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e(t)  \  x22(T)e(t)  =  a2p(t)  +  UQp(t)  +kZ=1ukp(t+Tk) 

B 


(77) 


and  this  equation  is  to  be  solved  subject  to  the  end- condition,  e(0)  =  e(T)  =  0, 
which  follows  from  Equation  (73).  Assuming  that  we  can  solve  this  equation 
for  e(.  ),  the  optimum  filter  impulse  response  for  the  combined  estimation  and 
detection  problems  is 


*  x5(T) 

h  (T-t)  =  e(t)  -  P(t) 

We  let 


(78) 


ai  =  -  x5(T)/x2(T) 


(79a) 


W2  -Ik*!  E 
a  ? 


x  2(T) 


(79b) 


It  is  clear  that  UU  is  a  positive  constant.  It  follows  from  (25)  and  the  definitions 
of  Xj-(T)  and  x^(T)  that  0^  must  be  positive.  Then  the  optimum  filter  impulse 
response  is 


h*(T-t)  =  e(t)  +  a^p(t)  a  >  0  (80a) 

where  e(t)  satisfies  the  equation 

2  n 

e(t)  -  0)  e (t)  =  a2p(t)  +  UQp(t)  +  kE  1  Ukp(t+Tk)  (80b) 

e  (0)  =  e(T)  =  0  (80c) 

There  remains  the  problem  of  selecting  the  multipliers  a^,  U  ,  •  ■  •  ,  U 

in  order  to  guarantee  that  the  specified  constraints  are  satisfied.  The  discussion 
of  this  problem  will  be  deferred  until  a  later  time.  The  main  point  to  make  here 
is  that  the  optimum  filter  is  made  up  of  parallel  processors:  in  one  path  is  a 
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matched  filter,  in  another  a  matched  filter  and  transversal  equalizer  are 
cascaded  with  some  sort  of  unstable  dynamical  system.  For  a  fixed  processing 
time  a  well-defined  solution  results,  but  for  the  case  in  which  T-*00,  an  unstable 
mode  arises  in  the  filter's  realization.  Furthermore,  we  really  are  interested 
in  obtaining  a  filter  impulse  response  on  [O,00)  and  this  is  not  possible  using 
the  state- space  technique  directly.  Presumably  we  could  use  the  method  out¬ 
lined  above  to  solve  the  problem  for  each  and  let  T^-*00  but  it  is  extremely 
difficult  to  study  the  behavior  of  the  optimum  solutions  as  a  function  of 
except  in  the  simplest  case  of  minimizing  the  detection  noise- to- signal  ratio. 
Therefore,  there  is  a  good  reason  to  re-examine  the  filter  design  problem  in 
the  frequency  domain  where  the  problems  of  stability  and  fixed  terminal  time 
do  not  arise.  We  do  not  mean  to  imply  that  the  optimal  control  formulation  has 
been  useless,  since  by  following  the  recipe  prescribed  by  the  minimum  principle 
we  have  succeeded  in  deducing  the  structure  of  the  optimum  filters.  This 
information  indicates  the  appropriate  inner  product  space  to  use  when  solving 
the  problem  in  the  frequency  domain. 
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IV.  Frequency  Domain  Solution  of  the  Design  Problem 


In  this  section  we  shall  use  the  L  (-00,  00  )  theory  of  Fourier  transforms  to 
transform  the  filter  design  problem  described  in  Equation  (40),  (4  1),  and  (42 ) 
into  the  frequency  domain.  Therefore  we  need  assume  that  h(.  ),  h  (.  ),  p(.  ), 
and  p(*  )  are  members  of  the  space 

o°  2 

L2 (-»,»)  =  {f(.  ):  J*|f(t)|dt<  » }  (81) 

_CO 

where  f(,)rnay  be  a  complex- valued  function.  Then  the  Fourier  Transform  of 
f(.  ),  is  well-defined  by 

F(uu)  =  J  f(t)e"jU)tdt  =Sf[f(t)]  (82a) 

_  00 

(82b) 

_  00 

where  the  M=M  signs  mean  ’’limit  in  the  mean”.  We  shall  have  occasion  to  use 
Parseval’s  Theorem 

OO  00 

J_oof(t)g(t)dt  =  F(u))G*(u))dm  (83) 

We  also  note  that  if  F(ud)  =^P[f(t)],  then  juuF(u))  =  P*[f(t)]  and  that  F  (uu)  =3*[f(-t)]. 

Using  these  relations,  the  filter  design  problem  stated  in  Equations  (4  0),  (41), 
and  (4  2)  becomes  the  problem  of  minimizing 


,2  J  w2|H(«u)|2 


dU) 


a 


“ET  : 


+  (1-a)- 


00  2  2 
[J  uu  H(U))P(uu)duu] 


1  H(ua)  l  2duu 

00 

[J  H(uu)P(uu)dw 


(84) 


0  =  a  =  1 


30 


subject  to  the  zero-bias  constraint 
00 

jxH(x)P(x)dx  =  0  (85) 

_  00 

and  sidelobe  constraints 


OO  00  j  (JjT  00 

-€kJ  H(ou)P(uu)du;  =  J  H(x)P(x)e  k  dx  =  j  H(x) P(x)dx  (86) 

.00  _  00  _  00 

k  =  l,  2 ,  .  .  .  ,  n 

The  inequality  constraints  are  as  troublesome  in  the  frequency  domain  as 
they  were  in  the  time  domain.  However,  we  can  use  the  same  trick  to  eliminate 
them  and  first  solve  a  problem  with  equality  constraints.  With  no  loss  of 
generality  we  can  set 


J  tu2H(uu)P(w)^  =  B2  (87) 

_  OO 

because  the  equations  are  homogeneous  in  H(x).  Furthermore  we  require  that 


j^T. 


J  H(u>)P(ou)e'  k^r  =  Yk  k=0,  1 _ _  n 


(88) 


n+i 

where  =  0  and  the  (n+1)- vector  y  belongs  to  the  S  c  E  where 


s  =  k'  -Vo  =  yk  =  Vo’  k=1’2 . n}  (89) 

If  we  consider  only  those  cases  for  which  a  >  0,  then  for  each  fixed  yES  we 
want  to  minimize 


— | - fj  UU  1  H(JL')  I  dx  + 

BE 


(1-0.)  B  E 


J  1  H(jl) |  dx] 


y  _  oo 

7  O 


subject  to  the  constraints 


00 

jxH(x)P(x)dx  =  0 

_  00 


(90) 


(91) 
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(92) 


j*  'JU2H('JU)P(uu)|^  =  BZ 
_  00 

J  HMPwg  =  yo 

_  00 

*  jujTkduu 

I  H(uu)P(U))e  =  yk  k=l,  2 . n 

_  00 

From  Parsevals  Theorem  we  know  that 


(93) 

(94) 


srJ..  lF(»)|2d»»j_.U(t)|2dt 


(95) 


Therefore  since  h(.  ).  h(*  ),  p(-  ),  p(*  )  belongs  to  L  (-<=, 00 )  then  H(uu),  uuH(uu)» 
P(uu),  a'P(uu)  belong  to  L  (-<»,«>).  Let  us  now  define  a  vector  spaceO^  over  the 
real  numbers  as 


=  [F(uu):  F(uu),  UUF(UU)  €  ]?(-»,«),  F(tu)P(U>)  =  F*(-0B)P*(-tt))}(96) 
The  condition  F(uu)P(uu)  =  F*(-UJ)P*(-uu)  is  equivalent  to  the  requirement  that  the 
complex  representation  of  the  compressed  pulse  be  real.  We  define  an 

inner  product  on  as  follows:  For  any  F(*  )  and  G(*  )e  , 

00  2 

<F(uu),  G(w)  >=  °  [uu2  +%^--  -5-J1L]  F(<u)G*  («>)**!"  (97) 

.00  a  y  2  IT 

'  o 

For  convenience  we  set 

,  2  _  (l-a)  b2E 

X  =  ~  “2”  (98) 

y  o 

and  note  that  this  is  a  positive  quantity.  It  is  easy  to  verify  that  Equation  (97) 
specifies  a  well-defined  inner  product.  Furthermore  it  can  be  verified  that 
the  pair  Ctf\>  <#  >  #>)  defines  a  Hilbert  Space.  We  can  then  imbed  the  minimiza¬ 
tion  problem  into  this  Hilbert  Space  and  minimize  Equation  (90)  which  becomes 
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a 


(99) 


-2-  IIhMI2 

BE 

2> 

where  |F(uu)||  =  <F(uu),F(uu)  >.  The  minimization  performed  subject  to  the  con¬ 

straints,  Equations  (91)  through  (94),  which  can  now  be  written  as 


<  H(uu),  >=  0  (10°) 

io  +X 

2  2 

<  H(u o),  — 2P*(U))>=B  (101) 

io  +X 


<  H(u»,  — i - zP*(uu)>  =  y  (102) 

uu  +X 

1  -j^T, 

<  H  (uu ) ,  — 2 — 2P*(UJ)e  >  =  yk  k=l,2,...,n  (103) 

uu  +X 


These  operations  are  well-defined  provided  the  functions  appearing  in  the  inner 
product  belong  to  the  space*^.  Since  H(uu)  and  P(uu)  belong  to^by  assumption, 
it  is  easy  to  show  that  the  functions  appearing  in  the  above  inner  products  belong 
to7f{  as  long  as  X  4  0.  Notice  that  X  =  0  if  and  only  if  <X  =  L  Referring  to 
Equation  (84),  this  case  corresponds  to  minimizing  the  estimate  variance  sub¬ 
ject  to  no  constraints  on  the  detection  signal- to- noise  ratio.  We  have  already 
given  physical  arguments  based  on  the  practical  problem  for  requiring  that 
a  ^  1.  Now  we  have  a  mathematical  requirement  that  this  condition  be  met. 

Next  we  define  a  set  L  to  be 


r  I00  UU 

L  =  1F(uj):  F(w)  =  a  -A, — ?  p*(u))  +  a  — — z  p*/,^ 

uu  +X  w  +X 


+  uo_T72  ***">  +  kE=luk_TT2ls'(",|e 

uu  +\  uu  +\ 


■j^T. 


(104) 
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where  a^,  a^,  UQ,  .  .  .  ,  Un  are  arbitrary  real  constants.  It  is  easy  to 
show  that  L  is  a  subspace  of  77[  .  Then  the  projection  theorem  for  functions 
in  a  Hilbert  Space  holds  [  14 ]  and  we  can  express  the  function  H(UU)€  y/la-s  follows 


H('JU)  =  Hl(uu  )  +  H±  (uo)  (105) 

where 

Ht  (uu  )€L  and 

<H^(uu),  F(u>)  >=  0  for  all  F(uu  )€L  (106) 

This  decomposition  of  H('JU)  is  unique.  Then  from  Equation  (99)  the  function  to 
be  minimized  is 


-2— II H  Ml  + 

BE  L 


a 

b2e 


IIh^MI2 


(107) 


All  of  the  functions  appearing  in  the  constraint  equations,  Equations  (100) 
through  (103),  belong  to  the  subspace  Land  are  therefore  orthogonal  to  (uu) . 
(Therefore  only  that  part  of  H('JU)  belonging  to  L  is  utilized  in  satisfying  the 
constraints.  )  Hence,  the  problem  reduces  to  minimizing  (107)  subject  to  the 
constraints. 


<  H  (uu),  -4£— , 2P*(uu)  >  =  0 

tu  +X 


<  H  (uu ),  g  2  >  =  B‘ 

uu  +\ 


<  T  2  )  >  "  yo 

uu  +X 


J_  .  -ja’Tk 

u u“+X‘ 


<  HL(m).  %  %P»(uu)e  ~  >  =  Yk  k=l,  2 


,  n 


(108) 

(109) 


(110) 

(HD 
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It  is  obvious  that  the  norm  can  be  made  smallest  by  choosing  H^(u,’)  =  0.  Then 

A 

the  optimum  solution  H(X')  belongs  to  L  and  therefore  is  of  the  form 


1  2  n  'JlTk 

H((U)  =  ~2 - j[a1jt«+a2UJ  4u0+l<F=lUke 

a)  +X 


(112) 


Since  the  problem  is  to  be  solved  for  a  fixed  a  and  y  with  0<  a  <  1,  v  2  0  and 
1  7  o  7  o 

2 

since  the  constants  a  a?  ,  U  ,  .  .  .  ,U  are  arbitrary  we  can  replace  a  by  X  a  , 

1  o  n  ii 

etc.  and  obtain 


,2  2  n  ”^aTk 

H(ul’)  =  ~2 - j  [a1j(iJ+a2a'  +u0+  kS=lUke  K]P>li('x’) 

uu  +X 


(113) 


Since 


X2  2  _  ,2  a2X 

~~2 — 2  aza)  “  a2  "  ~2 — 2 

JU  4*X  X  4*X 


(114) 


so  that  the  Fourier  Transform  of  the  optimum  filter  impulse  response  can  be 
written  as 

2  y  ^  ^  ~  j t_ 

1 1 (-r )  =  a2X  P*(oj)  +  -J— 2  [a  j:»+(u  -a2X  )+  .TjJ.e  ]  ^ ■  (  ) 

11  +X  (115) 


Since  a  a  u  ,  .  .  .  ,U  arbitrary  real  constants  we  can  relabel  them  at  our 
i^o  n 

convenience.  Similarly  the  constaint  points  can  be  relabelled  appropriately 
and  we  can  express  the  optimum  solution  as 


.2  n  -j  XT 

H(uu)  =  a  P*(uj)  +  — 2 - 2  [a2 +k=  n  Uke  (11&) 

a)  +X 

The  constants  a-,  a0,  u  ,  .  .  .  ,  U  are  to  be  chosen  so  that  the  solution  satisfies 
I  l  -  n  n 

the  constraint  equations  (108)  -  (111). 

Up  to  this  point  in  the  analysis,  the  function  H(JJ)  has  been  interpreted  as 

2 

the  -  Fourier  Transform  of  the  optimum  filter  impulse  response.  It  would 
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also  be  convenient  if  we  could  interpret  H  (uu)  as  the  transfer  function  of  the 
filter.  Let  us  define  the  functions 


f(t)  =  i  e_>  ' t' 


-00  <  t  <  CO 


n 


g(t)  =  ,z  u.  6(t-r,  ) 


(117a) 

(117b) 


Note  that  f(t)  can  be  interpreted  as  the  impulse  response  of  a  stable  but  un¬ 
realizable  filter,  while  g(t)  is  the  impulse  response  of  a  tapped  delay  line, 
(S(t)  represents  the  Dirac  delta  function).  Let  us  also  define  the  function 


q  (t)  =  aJp(t)  +  f(t)*[a2p(-t)  +  g(t)*p(-t)]  (118) 

where  ,f*,r  denotes  the  convolution  operation.  We  can  take  the  Fourier  Trans¬ 
form  of  q(t)  in  the  Generalized  Function  sense  to  obtain 

Q('ju)  =  a  P*((ju)  +  F(uu)[a2juu  +  G(u>)]  F%u)  (119) 

where  F(uu),  G(uu)  are  the  Fourier  Transform  of  f(t),  g(t).  In  other  words 

X2 

F  (lu) - 2 — t 

+X  (12  0a) 

n  -juur, 

G(uu)  =  k|_nuke  K  (120b) 

-  2 

But  the  function  H  (uo)  in  Equation  (115)  is  the  L  -  Fourier  Transform  of  h(t). 
Therefore  it  is  also  the  Fourier  Transform  K  (t)  in  the  Generalized  Function 
sense,  and  H(u u)  =  Q(uu),  in  other  words, 


H(uu)  =  a x  P^;(uu )  +  F(u))[a2ju)  +  G(i«)]P*(u))  (121) 

From  this  equation  we  can  draw  the  block  diagram  for  the  filter  as  shown  in 
Figure  2.  * 

*Note  that  in  Figure  2,  a^  =  0.  This  corresponds  to  a  result  shown  in  Chapter  V. 
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lu-M?444| 


Fig.  2.  Block  diagram  representation  of  the  optimum  mismatched  filter 


37 


Therefore  the  optimal  processor  for  estimating  the  target  range  and 
detecting  the  target  in  a  multiple  target  environment  consists  of  two  parallel 
processors,  one  of  which  is  simply  a  matched  filter,  the  other  a  matched  filter 
in  cascade  with  a  tapped  delay  line  and  a  bandlimiting  filter.  If  we  were  to 
study  simply  the  target  detection  case  in  clutter  we  would  have  chosen  the  cost 
weighting  factor,  a,  in  (84)  to  be  zero  and  we  would  not  have  needed  the 
normalization  (87)  but  rather  we  would  set  y^  =  1.  By  retracing  the  steps  used 
to  obtain  the  above  solution  (this  case  is  even  easier)  we  would  have  found  that 
the  optimum  filter  was  given  by  (121)  with  a^  =  0  and  F(d')  =  1.  Therefore  the 
filter  which  maximizes  the  detection  signal- to- noise  ratio  subject  to  sidelobe 
constraints  is  simply  a  matched  filter  followed  by  a  tapped  delay  line.  We 
have  therefore  established  the  optimality  of  the  classical  transversal  equalizer 
structure. 

In  other  cases,  in  which  a  ^  0,  the  cost  function  gives  some  emphasis  to 
the  estimation  performance  of  the  filter  and  the  optimum  structure  is  that 
shown  in  Figure  2.  Therefore  the  transversal  equalizer  by  itself  is  not  the 
optimum  filter,  when,  in  addition  to  detecting  the  presence  of  the  target  in 
clutter,  its  range  is  to  be  estimated  as  well. 

In  Section  VI  we  have  evaluated  the  performance  of  the  optimum  filter  for 
a  particular  transmitted  signal.  The  dynamics  of  the  various  filter  parameters 
are  shown  at  that  time.  For  our  purposes  at  this  time,  it  suffices  to  say  that 
the  solution  to  the  filtering  problem  has  been  found  using  Hilbert  Space  techniques 
in  the  frequency  domain.  All  components  in  the  filter  are  stable  but  the  band- 
limiting  filter  is  unrealizable.  No  questions  concerning  the  duration  of  the 
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processing  time  arise  in  this  case,  since  the  optimum  impulse  response, 

h(t)  =  J  l[H(uu)],  will  be  defined  on  the  whole  real  line,  -»  <  t  <  Therefore, 

the  frequency  domain  formulation  does  appear  to  have  advantages  over  the 

state- space  solution.  We  want  to  emphasize  again,  that  we  found  the  form  of 

the  solution  by  using  the  optimal  control  techniques  and  used  the  insights 

gained  from  that  solution  to  set  up  the  linear  space  and  the  inner  product 

which  led  to  the  mathematically  rigorous  result.  Furthermore  one  advantage 

of  the  state- space  approach  is  that  it  leads  directly  to  a  realizable  filter. 

The  problem  of  filter  synthesis  is  not  completely  solved,  however, 

since  we  have  yet  to  determine  the  constants  a_  ,  a~  ,  u  ,  .  .  .  ,  u.  .  In  the  next 
7  1  Z  -  n  n 

section  we  shall  show  how  these  parameters  can  be  found  by  solving  a 
quadratic  programming  problem. 
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V.  Reduction  of  the  Filter  Design  to  a  Nonlinear  Programming  Problem 

In  the  previous  section  we  showed  that  for  a  fixed  vector  y€S,  where 

s  =  {  YJ  _^kyo  ^  yk  ^  ^kyo’  k=±l,  ±2,  .  .  .  ,  ±n}  (122) 

the  function  which  minimized  the  combined  estimation/detection  cost  functional 


a 


B  E 

subject  to  the  constraints 


J  uj2|h(uj  )|2duj  +  |  H  (uj  )  |  2  doJ 


J  jujH(uj)P(u))|£=  0 


(123) 


(124a) 


J  uj2H(uu)P(uj  )^  =  B* 


(124b) 


3  j,JuTvi. 

H(uj)P(uj)e  -1—  =  k=0,±l,  .  .  .  ,±n 


was  of  the  form 


H(uj)  =  a  P*(<aj)  +  a2  T  P*(«>) 

uu  +X 


.2  n  -  JujTt_ 

+  — t - t  f  i  ?  ti.  e 


T— Z  L  k=-  nk 

uj  +X 


(124c) 


]P*(uj)  (125) 

The  constants  a  a  ,  .  .  .  ,  U  are  to  be  chosen  to  satisfy  constraint 

1  Z  -  n  n 


equation  (124).  In  this  section  we  propose  to  show  that  these  constants  can  be 
found  by  solving  a  nonlinear  programming  problem.  Before  we  do  this,  however, 
we  shall  show  that  the  (2n+l)-dimensional  constraint  set  can  be  reduced  by  a 
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factor  of  two  by  assuming  that  the  sidelobe  samples  are  taken  symmetrically 
about  Tq  =  0  and  that  the  sidelobe  constraint  function  is  symmetric.  Only  the 
latter  assumption  results  in  some  loss  in  generality  of  the  results.  First  we 
define  the  real-valued  functions 


X 


n 


)  SJUT, 


F r(X')  ”  al  +  2  "  2  ki-r^k00""'  k 

a  +X 

a^AX^  ^2  n 

Fi(lX))  =  "2  ~2  '  2  2  ki-  nUkSmXTk 

a  +X  A  +X 


and  note  then  that  Equation  (125)  can  be  written  as 


(126a) 

(126b) 


H('X)  =  [Fr(-ju)  +  jF  .(uu)]P*(u>) 
Now  we  make  the  assumption  that  the 
chosen  symmetrically,  i.  e. 


sidelobe  sampling  points  [t,  }n 

k  k=-n 


(127) 

be 


Tk  ”Tk 

This  assumption  results  in  no  loss  of  generality  of  the  method, 
us  to  rewrite  (126)  as 


Ft"”)  =  al  +  T— T  Iuo  +  kil  'V^k’00”^1 

(D  +A. 


a^'JUX  ‘ 


X 


Fi(X)  ‘  2  2  ’  2  2  k  =  l(Uk'U-k)sinX’Tk 

A  +X  cu  +X 

Let  us  define  new  variables  according  to 


(128) 


but  it  permits 


(12  9a) 

(12  9b) 


U  =  U  ,  U  =  a  ri  ,  =  a 

o  o  n+1  1  n+1  2 


Uk=Uk^-k’  ^V^k 


(130) 
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Then  the  function  F^(.  )  depends  only  upon  the  (n+2)- vector  u  while  F^(*  )  depends 
upon  only  the  (n+l)-vector  T)  .  These  vectors  are  to  be  chosen  to  satisfy  the 
constraint  equation,  (124).  From  (129)  it  is  obvious  that 


Fr(uu;u)  =  Fr(-uu;u)  (131a) 

F.(uu;ri)  =-F(-uu;p)  (131b) 

where  we  have  explicitly  denoted  the  respective  dependence  upon  uand  T|.  The 
optimum  solution  can  be  written  as 


H  (u?u ,  r\)  =  [Fr(lu;u)+jF.(UJ;i2)]P»!'(tJu)  (132) 

A 

We  now  substitute  for  H(.  )  in  the  constraint  equations  and  obtain  the  following 
set  of  constraint  equations  which  are  completely  equivalent  to  those  expressed 
in  (124): 


00  o 

r  uuF.(uu;r))|  P('ju)|21^  =  0 

loo  1  “  2ir 


(133a) 


J  uu2Fr(ui';u)|P(JU)|2g  =  B2 
—  00 

J  Fr(w;u)|P(w)|2^=yo 

-00  ^  " 

00 

[F  (uu;u  )cosJJT  -F.(uu;r))sinuuT,  ]|  P('ju)|  2^li  =  y,  (133d) 

^  r  k  i  K  Ztt  k 

—  oo 

These  equations  follow  from  the  fact  that  Fr(uu;\J)  and  |  P(tt))|  are  even  functions 

of  rr  (co  (t)  is  real)  while  F^(uu;^n)  is  an  odd  funct  ion  of  uu.  For  example  since 
2 

(JL'F  r  (au ;  u )  |  P(au)|  is  an  odd  function  of  <u 


(133b) 


(133c) 


42 


(134) 


and  the  zero-bias  constraint  is  reduced  to  Equation  (133a). 

The  cost  function,  (123)  depends  on  the  vectors  U  and  r\  according  to  the 
relation 


(135) 


.00  B  E  y 

7  o 

In  the  equations  for  Ff(*  )  and  (•  )  we  have  been  using  the  definition 


}  ^  s  (U) 


(136) 


a 


These  relations  have  been  derived  using  only  the  fact  that  T  ^  =  -T^.  It 
is  now  convenient  to  introduce  the  assumption  that  the  sidelobe  constraint 


function  is  symmetric.  This  condition  will  be  satisfied  for  the  so-called 


’’thumb- tack,f  constraint  function,  but  there  may  arise  cases  where  a  non- 
symmetric  sidelobe  structure  would  be  permissible.  We  would  handle  these 
cases  by  applying  the  constraints  symmetrically  which  would  implicitly  reduce 
the  degree  of  variation  in  the  multipliers  and  result  in  an  increased  cost.  We 
would  like  to  point  out  that  the  assumption  is  introduced  as  a  matter  of  con¬ 
venience  in  reducing  the  dimensionality  of  the  problem.  We  could  apply  the 
same  analysis  to  the  case  of  non- symmetric  constraint  functions,  and  thereby 
lose  no  generality  of  the  method.  The  symmetry  assumption  implies  that  the 
sidelobe  constraint  points  satisfy  the  condition 


yk  -  y_k  k=l,  2,  .  .  .  ,  n 


(13  7) 
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Using  this  condition  and  the  fact  that  r  k  =  -Tk>  it  is  easy  to  show  that  the  side- 
lobe  constraint  equation  (133d)  reduces  to  requiring  only  that 
00 

J  Fr(uj;u)cos'iJTk|  P(w)|  Z^1  =  yk  k=l,  2,  .  .  .  ,  n  (138) 

—  00 

Therefore,  for  each  (n+1)- vector  y  in  the  set  S,  where  now 


S  =  tx:  -€kyo  =  yk  =  6kV  k=l,  z,  .  .  .  ,  n}  (139) 

we  want  to  pick  the  constants  u,  T]  to  satisfy  the  reduced  set  of  constraint 
equations,  (133a),  (133b),  (I33c),  and  (139).  These  are  summarized  below  for 
conve  nience . 


J  uuF.(-juni)|P(uj|2|^  =  0 


(140a) 


J  a'2Fr(uj;u)|  P(uu)|2|^  =  B2 


(140b) 


f  F  (uu ;  u )  I  P(cju)|  2^£  =  y 


(140c) 


J  F^ (uj;u)cosa-'Tk|  P(oj)|  =  yk>  k=l,  2, 
_  oo  2tt 

and  the  associated  value  of  the  cost  functional  is, 


.  .  .  ,  n 

from  Equation  (135) 


(140d) 


J  (uu2+\2)[F  2(w;u)+F .2('Ju;ti)]|  P('JJ)|2duj  (141) 

B  E  -co 

Notice,  that  for  each  fixed  value  of  y€S,  that  only  the  zero-bias  constraint, 
Equation  (140a),  depends  upon  the  vector  r\.  All  of  the  other  constraints  are 
satisfied  by  the  appropriate  choice  of  u,  independent  of  T).  At  the  same  time, 
the  cost  function  is  increased  for  every  non-zero  value  of  r\.  From  (129b) 
and  (130), 
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(142) 


F.<u);t2)  = 


n 

E  TLsinUJT 


from  which  it  is  clear  that  zero-bias  constraint  can  be  satisfied  and  the  cost 


functional  minimized  by  choosing  T)  =  ^0,  i.  e.  F^(uu;p)  =  0.  Relating  this  con¬ 


dition  to  the  definition,  (130),  we  see  that  this  choice  of  n  is  indeed  possible 
by  picking 


(143) 


Therefore,  the  symmetry  of  the  sidelobe  constraint  function  permits  us  to 
reduce  the  dime nsionality  of  the  constraint  space  by  a  factor  of  2  and  we  no 
longer  need  concern  ourselves  with  the  zero-bias  constraint  since  this  will  be 
satisfied  automatically  by  the  above  choice  of  multipliers. 

Summarizing  these  results,  we  conclude  that  when  the  sidelobe  constraint 
function  is  symmetric,  the  optimum  solution  to  the  filtering  problem  is 


U(±)V)  =  F(uu;u  )P*(uu) 


(144) 


where 


£  n  u_  cosuur.  +u  , , 
c  =  0  k  k  n+1 


(145) 


and  the  performance  of  this  filter  is  measured  by  the  cost  functional 


00 


(14  6) 


where 


2  _  (1-ct)  BZE 
a  ‘  2 


(147) 


and  where  the  (n+2)- vector  u  is  chosen  to  satisfy  the  constraints 
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(148a) 


00 


J  'JU2F(uj;u)|  P(uu)|  =  b2 


00 


I  F(W;u)|P(uJ)|2g.  yo 

-  00 


(148b) 


00 


(148c) 


_  00 


In  a  crude  sense,  we  could  imagine  picking  a  vector  in  the  set  S,  solving 
the  preceeding  (n+2)  constraint  equations  for  the  (n+2)  vector  u  which  can  then 
be  used  to  obtain  a  value  of  the  cost  function  for  that  particular  vector. 

This  could  be  done  for  all  possible  y-vectors  in  S  and  the  particular  y  which 
minimized  the  cost  could  be  found.  In  the  next  few  paragraphs  we  shall  show 
how  this  search  procedure  can  be  performed  in  a  systematic  way  using  quad¬ 
ratic  programming  techniques.  We  begin  the  analysis  by  first  manipulating 

the  constraint  equation  (148).  From  (147)  we  notice  that  for  each  fixed  value 
2 

of  y  ,  X  is  a  fixed  positive  number.  Therefore,  if  Equation  (148b)  holds, 
the  following  equation  holds. 


7 


oo 


(149) 


Then  if  (148b)  holds,  constraint  (148a)  is  equivalent  to  requiring  that 


(150) 


Therefore  for  each  ^r6S  we  solve  for  the  vector  u  which  satisfies  the  constraints, 
Equations  (148b),  (148c),  and  (150),  namely 


00 


(151a) 


-  00 
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00 


J  F(a);u)cosUL'T  I  P(uu)|z|ii!  =  y  k=l,2,...,n 


2ir 


(151b) 


2.,2_  ,  - 1  2dau  _  r2[l|(l-a)  Ej 


J  (tu  +\  )F  (uu;u)|  P(x)| 


a  y. 


(151c) 


and  the  associated  cost  will  be 


a 


o  o  n 

J  (u  +\  )F  (uu;u)|  P(uu)|  ^duu 
B“E  -o° 


where 


(152) 


F(«,;u)  =  ~Z~2  Ao  ukcos;,;Tk  +  u„+l 

0)  +X 


(153) 


Now  define  the  complex  vector- valued  function  V(jo)  as 


v(»)  =  ^)*x- 


-j'd)+x 


1 

CO  SOOT  ^ 
C  O  SUUT  - 


CO  SUIT 

n 

2  2 
U)  +X 


\ 

and  the  (n+2)  x  (n+2)  non- negative  definite  matrix 


(154) 


(155) 


Here  the  asterisk  denotes  complex  conjugate  while  the  prime  denotes  vector 
transpose.  It  is  easy  to  show  by  straight  forward  manipulations  that  the  con¬ 
straint  equations,  Equation  (151)  can  be  written  more  compactly  as 
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ru  = 


yo 

yl 

y2 


y 


a 


B  [i+'iSL.  s.) 
a  v 


(156) 


and  for  each  value  of  u  satisfying  these  constraints  the  associated  cost  is 


f(u)  =  ^V(u_'ru) 

B  E 

Since  the  vector  then 


(157) 


y  =  y,  =  y  k=l,  2,  .  .  .  ,  n 


"k7  o  yk  k7o 

Because  the  constraints  on  U  are  actually  these  inequality  constraints  there 
is  a  certain  amount  of  freedom  in  the  selection  of  the  u-vector,  We  would 


(158) 


like  to  find  the  particular  value  of  u  within  prescribed  class  which  minimizes 

2 

the  cost  function.  Notice  that  for  each  value  of  y  ,  X  is  fixed  number  and 

o 

therefore  the  vector  V('JU)  is  a  completely  known  function.  Consequently  the 
matrix  T  is  fixed  also,  but  it  changes  as  y^  changes.  Let  us  denote  this 
dependence  explicitly  by  writing  V(uu;yo)  and  r(yQ).  Then  for  each  value  of  y^ 


we  want  to  find  the  vector  u  which  minimizes  the  cost  function 

oa2(y) 

f(u;y  )  = - £ -  •  (U  r (y  )u ) 

°  BE  “  °“ 


(159) 


subject  to  the  constraints 
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■6lyo 


-6 


2yo 


-6  y 
rr  o 


7 - [1 +Q£1  E] 

X  W  a  y 


=  Hy0)u  = 


Vo 

Vo 


B2 


n  o 

2  rn(1-a)  E  l 
TT. — ;ll+“a—  ^TJ 


X“(y0) 


(160) 


which  is  a  nonlinear  programming  problem.  Since  (160)  repre  se  nts  linear  in¬ 
equality  constraints,  Hildreth's  Method  [9]  can  be  used  to  transform  this 
optimization  problem  to  a  quadratic  programming  problem  of  the  form:  Find 
the  vector  x  =  0  which  minimizes  the  cost  function  c  'x  +  x '  Dx.  Efficient 
computational  algorithms  which  guarantee  convergence  to  an  optimum  x  in  a 
finite  number  of  steps  are  available  to  solve  this  latter  problem  flO]. 

Therefore  for  each  y  ,  we  can  solve  for  the  minimizing  u  =  U  (y  )  which 
has  associated  with  it  the  cost  f(u(yQ);yo).  We  then  use  conventional  search 
techniques  [ll]  to  search  over  all  possible  values  of  y^  to  find  the  pair  (u*  ,  y  *) 
where  V*  =  U(Y0*)  which  minimized  f(u(y^),yo)  treated  as  a  function  of  y  . 

This  technique  can  be  implimented  on  the  computer  in  a  relatively  straight¬ 
forward  way.  The  major  difficulty  lies  in  the  evaluation  of  the  matrix  F(y^) 
defined  in  Equations  (154)  and  (155).  For  reasonable  choices  of  the  transmitted 
signal,  the  evaluation  can  often  times  be  performed  analytically  without  recourse 
to  the  digital  computer. 

n+1 

Therefore  the  optimum  set  of  multipliers  {u^}  and  the  filter  pole 
location,  described  by 
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(161) 


2  _  (1-a)  B2E 
~  ~ a  — T 

y  „ 


can  be  found  for  any  weighting  factor  a.  The  resulting  optimum  filter  is 


(162) 


H((ju)  =  F(u>;u)P*(U>) 


where 


u,  cos'JJT,  +  u  , , 
:=0  k  k  n+1 


(163) 


The  actual  estimate  variance  and  detection  signal- to- noise  ratio  can  then  be 
evaluated  for  any  value  of  the  weighting  parameter  Ct.  We  would  construct 
a  table  of  such  values.  Then  for  a  prescribed  application  the  desired  detection 
signal- to- noise  performance  would  be  specified  and  we  would  reference  the 
table  we  had  constructed  to  find  the  appropriate  value  of  a.  With  this  weighting 
we  could  then  solve  for  the  correct  values  of  the  multipliers  and  thereby 
generate  the  filter  which  gives  the  desired  detection  performance  and  yields 
the  smallest  estimation  variance.  In  the  next  section  we  shall  apply  this 
methodology  to  a  practical  design  problem. 
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VI.  Numerical  Results 


In  order  to  perform  a  numerical  study  of  the  dynamics  of  the  parameters 
involved  in  the  optimum  mismatched  filter,  it  is  necessary  to  evaluate  the 
matr  ix 


J  V*(w;yo)V>;yo)g 


(164) 


where 


V(:u;yo)  = 


P(x)\2 

-ji+X 


1 

c  o  suur  i 
co  suT  ^ 


COSiT 

r 

2  2 

(D  +X 

“2 — 


X 


(165) 


,  Z  _  (1-0.)  b2e 

K  ~  ~d  T  (166) 

y 

7  o 

Conceptually  this  evaluation  can  always  be  done  since  once  the  transmitted 
signal  is  specified,  P(d)  could  be  found  and  the  integrations  performed.  This 
could  be  a  costly  exercise  if  recourse  to  a  digital  computer  need  be  made, 
since  the  matrix  is  of  large  dimension  and  would  have  to  be  evaluated  for 
possibly  many  values  of  y  .  For  our  purposes,  we  mainly  want  to  demonstrate 
the  feasibility  of  the  design  method  and  we  shall  therefore  restrict  our  interest 
to  a  class  of  signals  for  which  the  matrix  computation  can  be  done  analytically. 
We  motivate  our  choice  by  noting  that  the  vector  V(x;y  )  could  be  written  more 
generally  as 
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1 


,r/  ,  P(s)X 

^(s;y0)  =Ti+r 


i  sTi,  i  -STi 
e  +  ^ 


2 

h  ~+ i? 


1  _  ST2  .  1  -ST, 


-ST  ,  -  ST 

1  n  ,  1  n 

2®  +  7  6 

(-  s+X  )(s+X ) 

.  1T~  . 


(167) 


where  we  have  merely  replaced  juu  by  the  complex  variable  s  =  a+juu.  Then 
the  integral,  Equation  (164)  can  be  viewed  as  the  value  of  the  contour  integral 
of  the  complex  function 


r(yG)  • 


l  rJ“ 

72  J  x(- 

X  -j00 


s;y 


)V'( 


O  — 


s;y0) 


ds 


(168) 


where  the  contour  is  chosen  in  the  s-plane  along  the  juu  axis.  Then  if  we  choose 
our  time  functions  so  that 

m  a 


P(s)  =  E 


m 


m=l  (s  +  s  )(s  +  s  -) 
m  m  7 


(169) 


then  we  can  actually  evaluate  the  integral  using  the  residue  theorem.  Our 
results  are  further  restricted  to  the  particularly  simple  signal  for  which 


(s  +  s  )(s  +  s 
V  P  P 

We  shall  choose  the  constant, 

i.  e.  ^ 

E=J  |p(ja,)|2|^  = 
_  00 


a, 


1 


(170) 

to  normalize  the  transmitted  signal  energy, 

(171) 
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and  we  can  then  compute 


B2  =  J  u>2|P(jtu)i2^  (172) 

_  00 

which  now  represents  the  second  moment  bandwidth  of  the  signal.  We  can 

therefore  change  the  signal  bandwidth  by  varying  the  pole  location  s  .  We 

P 

then  used  (170)  in  (168)  and  applied  the  residue  theorem  as  proposed.  This 

left  us  with  a  closed  form  expression  for  TT(y^)  in  terms  of  y^,  and  this  was  easily 

programmed  on  the  digital  computer. 

If  we  represent  s  in  (170)  by  CF  +  j(JU  ,  the  compressed  pulse  for  a  matched 
P  P  P 

filter  is  given  by 


-at  a 

co(t)  =  e  P  (cos uu  t  +  — —  sin  X  t),  t  ^  0  (173) 

p  ju  p 

P 

withcp(-t)  =  CD (t ) .  The  mean-square  bandwidth  of  the  signal  is  easily  calculated 
to  be 


2  Z  Z 
B  =o)  4*  a 

P  P 

By  diffe r e ntiating  (173), 
found  to  be  located  at 


(174) 

the  subsidiary  peaks  of  the  compressed  pulse  are 

=  k7r/(JL  ,  |k|  =  1,2 . with  the  peak  magnitudes 

tra  P 


given  by  |cD(t,^)|  =  exp(-)kl  —  — ) .  The  choice  of  a  =  471^  results  in  four 

' P  A  - 1  P  P 

sidelobes  within  a  single  time-constant  (  =  a  ),  the  first  peak  magnitude 

P 

being  -  2.3Dbwith  respect  to  the  main  peak.  The  compressed  pulse  and  its 

spectrum  for  a  signal  a  =  1,  X'  =  4t t  are  shown  in  Figure  3. 

P  P 

To  illustrate  the  techniques  developed  in  this  paper,  constraint  levels 
corresponding  to  -10  Db  and  -20  Db  sidelobes  were  considered.  In  each  case, 
the  constraint  wa  s  applied  at  51  evenly  spaced  points  in  the  c losed  inte rval  [0.  25 , 2 ,  0]. 
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Fig.  3a.  Matched  filter  compressed  pulse. 
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Fig.  3b.  Compressed  pulse  spectrum. 
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Relating  this  problem  to  the  development  in  Section  5,  we  have  chosen 

T  =  +  T,,  6  =  6,  n=l,  .  .  .  51=N  where  T  =  0.  25  and  T  =2.0  . 

n  N-l  N  1  1  n  1  N 

For  a  fixed  a-value,  the  quadratic  programming  problem  determined  by  (159) 
and  (160)  is  solved  for  several  values  of  y^.  The  minimum  cost  function  is 
plotted  vs  y^  in  Figure  4  for  a  =  0.  5  .  For  any  arbitrary  a,  the  minimizing 

value  of  y^  can  be  determined  by  such  graphical  means  and  will  be  denoted 
by  y  *(a).  The  solution  of  the  quadratic  programming  problem  specified  by 
OL.y  *(CL)  will  be  denoted  by  U*(a).  The  bandwidth  parameter  for  this  optimum 
solution  can  be  calculated  from  (161)  with  y^  =  yQ*(&)  and  will  be  denoted  by 
X*(a) 


Recall  that  the  cost  function  given  by  (159)  is  a  weighted  sum  of  the 

estimation  and  detection  performance  terms.  We  wish  to  determine  the 

performance  values  separately  for  the  optimum  processors  H*((JU  ;a  )  which 

are  obtained  from  ( 16  2  )  and  (163)  with  U  =  u*(a)  and  X  =  X  *  (& ) .  The  resulting 

^2 

estimation  performance,  £  -  — j -  ,  is  obtained  as  a  function  of  a  and  is 

a  mf 

calculated  by  the  formula 


£  (CL)  =  B 


-|_00 

J  uj  [F*(uo;a)]2 


(175) 


where 


N 


F*(uu;a)  = 


a j 


-x — Uf:cosa'T.  +u*Tl. 
2  k=o  k  k  N+l 


(X*(a))+1  ^ 

The  detection  performance,  j£)  -  ~ p^£  >  can  he  calculated  by 


(176) 


3  (a)  = 


E 

[y  *(a)]Z 


J+  [F*(u»] 


_  CO 


2  d  au 


(177) 
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F  rom  (175)  and  (177),  the  optimum  trade-off  between  detection  and  estimation 
performance  can  be  determined  by  constructing  the  graph  given  by  the  set  of 
ordered  pairs  {  (  $  \a),  £  \a)):  0<a<ll  .  For  the  particular  signal  chosen, 

two  such  graphs  c  or  re  sponding  to  two  sidelobe  constraint  levels  are  shown 

-i 

in  Figure  5.  For  any  given  level  of  detection  performance  (^corresponds  to 
effective  SNR  loss  due  to  mismatching),  the  graphs  of  Figure  5  indicate  the 
best  possible  range  estimation  which  can  be  attained  subject  to  the  appropriate 
sidelobe  constraint  (and,  of  course,  the  original  assumptions  on  receiver 
structure). 

For  purposes  of  comparison,  the  performance  of  a  well-known  sub-optimal 
processor  using  the  c lutter- r eje ction  filter  was  calculated.  The  clutter  re¬ 
jection  filter  is  specified  by  [12],  1 13 ] 


H  (x)  - 

cr  D  +  |  P(a')| 


(178) 


The  sidelobe  level  achieved  by  this  filter  depends  on  p.  For  our  particular 
choice  of  signal,  the  compressed  pulse  and  the  performance  of  H  ^  can  be 
analytically  determined  as  a  function  of  p  using  standard  La  Place  transform 
techniques.  The  analysis  is  somewhat  lengthy,  however  so  it  is  omitted. 

The  performance  of  the  clutter  rejection  filter  for  several  sidelobe  levels  is 
indicated  on  Figure  5  for  comparison  with  the  optimal  values. 

A  typical  compressed  pulse  and  its  spectrum  resulting  from  the  optimum 
filter  for  a  =  0.  5  is  shown  in  Figure  6.  The  sidelobe  constraint  level  was 
-20Db,  cor  re  sponding  to  a  peak  sidelobe  magnitude  of  0.  1  .  The  sidelobe 

structure  of  the  compressed  pulse  actually  exceeds  this  level  by  a  slight 
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COST  FUNCTION 


Fig.  5.  Detection  and  estimation 
performance  of  optimum  mis¬ 
matched  filter. 


CLUTTER  REJECTION  FILTER 
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NUMERICAL  LABEL  “20 
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IN  dB 


0  5  10 

DETECTION  LOSS  (dB) 
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PULSE  SPECTRUM  (dB) 


Fig.  6a.  Optimum  mismatched 
filter  compressed  pulse. 


Fig.  6b.  Compressed  pulse 
spectrum. 
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amount.  This  is  due  to  the  fact  that  the  constraint  is  imposed  at  only  a  finite 
number  of  points.  In  between  these  constraint  points,  the  sidelobe  level  may 
exceed  the  desired  value.  This  effect,  if  observed,  may  be  reduced  by  simply 
taking  more  constraint  points  in  the  violated  region.  This  procedure  does 
not  necessarily  increase  the  number  of  required  taps  in  the  delay  line  realization 
of  the  optimum  filter  as  may  appear  to  be  the  case.  A  property  of  quadratic 
programming  problems,  when  interpreted  in  the  context  at  hand,  states  that 
the  gain  of  a  particular  tap  will  be  zero  unless  the  constraint  cor  re  sponding 
to  that  tap  (Equation  ( 16 0 ) )  is  satisfied  with  equality.  In  other  words,  if 
U*  ,  the  gain  of  the  tap  corresponding  to  a  second  delay  (advance)  is  non¬ 
zero,  then  the  value  of  the  compressed  pulse  at  seconds  must  be  either 

±  €  y  * .  In  most  cases,  the  number  of  times  this  condition  occurs  is  deter- 
n'  o 

mined  more  by  the  number  of  subsidiary  peaks  in  the  compressed  pulse 
structure  than  by  the  number  of  constraint  points  chosen.  For  the  case 
studied  here,  the  number  of  non- zero  tap  gains  (single- sided)  varied  from 
one  to  four. 
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VII.  Conclusions 


The  results  of  the  preceding  section  show  that  the  techniques  developed 
in  this  paper  can  produce  a  significant  level  of  sidelobe  reduction  whi  le  mini- 
mizi  ng  loss  in  effective  signal- to- noise  ratio  and  range  estimation  accuracy. 
Specifically,  we  have  considered  the  problem  of  minimizing  a  weighted  sum 
of  the  estimation  and  detection  performance  over  a  broad  class  of  receivers 
subject  to  a  sidelobe  constraint.  We  have  obtained  the  solution  to  this  pro¬ 
blem  for  an  extremely  large  class  of  signals,  specifically  those  with  finite 
energy  and  finite  mean- square  bandwidth  and  with  a  symmetric  spectrum; 
furthermore,  the  structure  of  the  optimum  filter  is  quite  simple. 

For  the  case  in  which  detection  performance  alone  is  considered,  the 
optimal  filter  consists  of  a  matched  filter  followed  by  a  transversal  filter, 
the  so-called  transversal  equalizer.  We  have  therefore  established  the 
optimality  of  this  well-known  receiver  structure  for  the  case  of  target  detec¬ 
tion.  In  the  more  general  case  where  estimation  accuracy  is  considered,  the 
transversal  equalizer  is  no  longer  the  optimum  processor;  further  filtering 
is  required.  The  nature  of  this  additional  filtering  is  such  that  the  width  of 
the  main  peak  of  the  compressed  pulse  is  increased  as  emphasis  is  shifted 
to  estimation  performance,  a  somewhat  surprising  result  since  widening  of 
the  central  peak  is  often  considered  synonomous  with  a  reduction  in  estimation 
accuracy.  In  Figure  7,  the  compressed  pulses  of  the  optimum  filters  for 
three  different  values  of  or  are  shown.  As  the  estimation  performance  is 
emphasized  (a  is  increased),  a  widening  of  the  pulse  is  observed.  In  order 
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Fig.  7.  Comparison  of  compressed  pulse  main  lobe. 
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to  understand  this  phenomenon,  the  basic  equation  for  the  normalized 

a2 


estimation  variance, 


mmf 


orrf 


is  written  in  the  following  form 


T  T2 
,  mmf  mmf 


T2  IT2  ' 
mmf  mf 


mmf 

Pmf 


T  ^  -  l/4^  —  y(0)/y(0)  where  y(t)  is  the  compressed  pulsed  of  a  mismatched 

filter,  h(t).  0.'  r  is  the  mean- square  bandwidth  of  h(t).  Note  that  T  r  is 

mmf  mmf 

a  measure  of  the  width  of  the  compressed  pulse.  For  a  matched  filter  this 
pulsewidth  is  written  as  and  is  given  by 


— z* 

t  f  ■  w\T  -z  where  a rn^s  the  mean- square  bandwidth  of  the  signal 
mf 

(and  the  matched  filter).  As  the  weighting  is  shifted  from  detection  to  estima- 

P 


tion  performance,  the  detection  mismatch  loss, 


mmf 

mf 


,  naturally  increases. 


Since  the  pulse  width  ratio,  T  -/T  r  also  increases  with  n  as  indicated 
r  mmf  mf 

by  Figure  7,  the  only  way  the  estimate  variance  can  possibly  decrease  is  for 
the  filter  bandwidth,  a  r,  to  be  reduced.  The  necessary  bandwidth  reduction 

is  accomplished  by  decreasing  the  parameter  \  *.  The  table  in  Figure  7  shows 
for  the  three  pertinent  values  of  a.  The  characteristic  of  the  optimum  filter, 
at  least  for  the  signal  considered,  we  note,  is  that  good  estimation  performance 
is  achieved  by  allowing  the  compressed  pulse  to  spread  by  simultaneously 
reducing  filter  bandwidth.  In  sharp  contrast,  the  clutte  r- rejection  filter 
behaves  in  the  opposite  manner.  The  compressed  pulse  of  this  filter  is 
actually  narrower  than  the  matched  filter  output,  but  the  bandwidth  is 
substantially  greater,  [  13 ] - 
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The  most  restrictive  assumption  made  in  this  paper  was  that  the  com¬ 
pressed  pulse  (complex  representation)  be  real.  This  restriction  excludes 
the  possibility  of  a  non- symmetric  signal  spectrum.  Also,  the  possibility 
of  treating  perturbations  caused  by  an  imperfectly  constructed  matched  filter 
or  by  a  distorted  input  signal  is  denied  since  both  these  situations  lead  to  a 
complex  compressed  pulse  in  general.  If  the  compressed  pulse  is  allowed 
to  be  complex,  the  performance  criteria  (estimate  variance  and  detection 
mismatch  loss)  must  be  modified  to  include  a  frequency  bias  term.  The 
more  general  formulas  may  be  found  in  [3].  Furthermore,  the  sidelobe 
constraints  will  be  quadratic  in  nature  rather  than  linear.  An  equivalent 
problem  for  this  more  general  case  would  probably  have  to  involve  a  more 
general  convex  programming  algorithm  and  an  additional  search  over  the 
frequency  bias  term. 
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13.  ABSTRACT 

In  a  multiple  target  environment  a  radar  signal  processor  often  uses  weighting  filters  which  are  not 
necessarily  matched  to  the  transmitted  waveform.  In  this  paper  expressions  for  the  mean-square  range- 
estimation  error,  the  detection  signal-to-noise  (SNR)  and  the  effects  of  sidelobes  are  derived  in  terms  of  the 
impulse  response  of  an  arbitrary  mismatched  filter.  It  is  desired  to  find  that  impulse  response  which  results 
in  the  minimum  range  estimate  variance  subject  to  preassigned  constraints  on  the  sidelobes  and  the  detection 
SNR.  This  optimization  problem  is  first  formulated  in  state-space  in  which  the  optimal  control  law  is  sought. 
Pontryagin’s  maximum  principle  is  used  to  obtain  necessary  conditions  for  the  optimum  impulse  response, 
from  which  it  is  possible  to  deduce  the  structure  of  the  optimum  filter.  Certain  mathematical  details  which 
detract  from  the  rigor  of  the  time  domain  formulation  are  resolved  by  formulating  the  problem  in  the  frequency 
domain  and  applying  Hilbert  Space  techniques.  It  is  shown  that  for  the  problem  of  detecting  the  radar  target 
and  estimating  its  range,  the  optimum  filter  is  a  modified  transversal  equalizer.  If  only  the  detection  function 
is  to  be  performed  the  optimum  filter  reduces  to  the  transversal  equalizer.  This  establishes  the  optimality 
of  this  important  practical  device  as  the  solution  to  the  radar  detection  problem  in  a  multiple  target  environ¬ 
ment.  The  tap  weights  and  spaces  of  the  delay  line  as  well  as  certain  other  parameters  upon  which  the  solution 
depends  can  be  found  by  solving  a  non-linear  programming  problem.  Numerical  results  are  given  for  an  inter¬ 
esting  class  of  transmitted  waveforms  which  shows  the  tradeoffs  of  the  various  filter  parameters. 
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